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Abstract We address the question of understanding the effect of the underlying network topology on 
the spread of a virus and the dissemination of information when users are mobile performing independent 
i random walks on a graph. To this end we propose a simple model of infection that enables to study 

the coincidence time of two random walkers on an arbitrary graph. By studying the coincidence time of 
a susceptible and an infected individual both moving in the graph we obtain estimates of the infection 
' probability. 

p ^ , The main result of this paper is to pinpoint the impact of the network topology on the infection 

probability. More precisely, we prove that for homogeneous graph including regular graphs and the classical 
Erdos-Renyi model, the coincidence time is inversely proportional to the number of nodes in the graph. 
' We then study the model on power-law graphs, that exhibit heterogeneous connectivity patterns, and show 

the existence of a phase transition for the coincidence time depending on the parameter of the power-law 
of the degree distribution. 

QQ ■ 1 Introduction 

In recent years there have been a surge of hand-held wireless computing devices such as PDAs together 
CO ' with the proliferation of new services. These portable computing devices are equipped with a short-range 

i wireless technology such as WiFi or Bluetooth. Despite providing a great deal of flexibility this ability 

to wirelessly connect to other devices, and to transfer data on the move, attracted the attention of virus 
. writers who exploit such features for lauching computer-virus outbreaks that take advantage of human 

mobility piTT]. 

Over the past couple years, there have been indeed reports of malicious code that take advantage of 
KJi \ bluetooth vulnerabilities such as the Cabir worm that was detected during the World Athletics Champi- 

i onship [12] and another at a company that has been reported by Comm Warrior [ISj. Despite their small 

\ scales, these incidents bode more threats taking advantage of events and locations where individuals gather 

in close proximity [19l [18] . 

In much of the literature on mathematical epidemiology, the members of the population are assumed 
to occupy fixed locations and the probability of infection passing between a pair of them in a fixed time 
interval is taken to be some function of the distance between them. Mean-field models are a special 
case in which this function is a constant [S]. In this work, we consider a different model in which the 
agents are mobile and can only infect each other if they are in sufficiently close proximity. The model is 
motivated both by certain kinds of biological epidemics, whose transmission may be dominated by sites at 
which individuals gather in close proximity (e.g. workplaces or public transport for a disease like SARS, 
cattle markets for foot-and-mouth disease, etc.) and by malware spreading between wireless devices via 
Bluetooth connections, for example. 

Related work. Here we briefly describe some of the relevant related work on modelling epidemic spreading 
in mobile environments. To our knowledge the first attempts to model virus spreading in mobile networks 
relies on the use of a non-rigorous mean-field approximations (similar to the classical Kcphart- White 
model [To]) that incorporates the mobility patterns of users. In [16], the authors derive a threshold for 
the persistence of the epidemic by computing the average number of neighbours of a given node. Using a 
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similar approach but with different mobihty patterns, Nekovee et al. [HI [T5] explore the evolution of the 
number of devices that are infected in terms of the contact rate betwen users. 

A related line of work studying the dissemination of information in opportunistic networks [3J focuses 
on the following analogous problem: Suppose that all individuals are interested in a piece of information 
that is initially held by one user. The information is transmitted between users who happen to be close to 
each other. As in the case of static networks [17 , one may be interested in the time it takes for the rumour 
to be known to all users. To this end we need to understand how information is transmitted between an 
informed and an ignorant user. Our work gives some insight on the impact of the network structure on 
the likelihood of successfully transmitting the rumour. 

Our contribution. In contrast to the previous work which has focused on Euclidean models and homo- 
geneous mobility patterns, in this work we consider a model wherein the different locations that a user 
can reach have varying popularity. 

More precisely, we consider a simple and stylised mathematical model of the spread of infection as 
follows. There is a finite, connected, undirected graph G = {V, E) on which the individuals perform 
independent random walks: they stay at which vertex for an exponentially distributed time with unit 
mean, and then move to a neighbour of that vertex chosen uniformly at random. The infection can pass 
from an infected to a susceptible individual only if they are both at the same vertex, and the probability 
of its being passed over a time interval of length r is 1 — exp(— /3t), where /3 > is a parameter called 
the infection rate. We shall consider a single infected and a single susceptible individual and ask what 
the probability is that the susceptible individual becomes infected by time t. This probability has been 
studied in the case of a complete graph in [6]. Here, we extend their results to a much wider class of 
graphs. 

It is simplistic to consider just a single infective and a single susceptible individual. Nevertheless, 
insights gained from this setting are relevant in the "sparse" case, where the number of both infected 
and susceptible individuals is small and inter-contact times are fairly large. In that case, it is not a bad 
approximation to consider each pair of individuals in isolation. The "dense" setting will require quite 
different techniques and is not treated here. 

The rest of the paper is organised as follows. In Section [21 we present our model and the family of 
networks we will consider. Besides we state our main results that relate the coincidence time of the two 
walkers to the stationary distribution of a random walk on a graph. In Section [3] we give a detailed proof 
of our main result on the probability of infection for regular graphs, the Erdos-Renyi graph and power-law 
networks. Section |4] summarises our contribution and suggests further extensions of our work. 



2 Models and results 

We now describe the model precisely. Let Xt^Yt G V denote the positions of the susceptible and infected 
individuals respectively at time t. We model {Xt,t > 0) and {Yt,t > 0) as independent continuous-time 
Markov chains (CTMCs) on the finite state space V, with the same transition rate matrix given by 

[ degrL(.) if y) e E, 

Qxy = I \iy ^ X and (x, y) ^ E, 

[ — 1 if y = cc 

where degree(x) is the number of neighbours of x (nodes y such that [x, y) G E) in the graph G. 
We define the coincidence time up to time i, denoted T(i), as the total time up to t during which both 
walkers are at the same vertex, i.e., 

T{t) = [ l{Xs - Ys)ds. (1) 
Jo 

Let 7(t) denote the probability that the initial susceptible becomes infected by time t. Then, conditional 
on r(i), we have 

7(t) = 1 - exp(-/3r(t)), (2) 

where /3 > is the infection rate. 
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We are interested in estimating the coincidence time r(i) and the infection probabihty ^(t) for different 
families of graphs. 

Observe that the Markov chains Xt, Yt have invariant distribution tt given by 

^ ^ degree(a;) 

E^GV degree(u) 

and that they are reversible, i.e., n^q-xy — i^y^yx for all x,y . 

We consider the case when these chains are started independently in the stationary distribution and 
provide estimates on the coincidence time and the infection probability, for arbitrary graphs. 

Theorem 1. Suppose Xq and Yo are chosen independently according to the invariant distribution tt. Then, 
we have 

E[T(i)] = ^[^(^)] ^ 1 - exp(-/3i ^l) ■ 

Proof. Observe that, for all s > 0, 

nxs = ys) = Y. = Ys = v) = Y 

because Xs and Ys are independent, and are in stationarity. Hence, it is immediate from ([T]) that 
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Next, taking expectations in ([2]) with respect to the conditioning random variable T{t), we have 

E[7(t)] = 1 - E[exp(-/3T(t))] 
< 1 - exp(-/3E[r(i)]) 

= l-exp(^-/3iE^^ 

where the inequality follows from Jensen's inequality. □ 
We now introduce some terminology and define some examples of graph models that we shall consider. 
For two functions /(•) and g{-) on the natural numbers, we write f(n) ~ g{n) to mean that their 
ratio tends to 1 as n tends to infinity. We write /(n) = 0{g{n)) if f(n)/g{n) remains bounded by a finite 
constant, /(n) = o{g{n)) if f{n)/g{n) tends to zero, and f{n) = n{g{n)) ii g{n) = 0{f{n)). For a sequence 
of events An indexed by n G N, we say that they occur with high probability (whp) if P(j4„) tends to 1 as 
n tends to infinity. 

Examples 

Complete graphs. Consider the complete graph on n nodes, namely the graph in which there is an edge 
between every pair of nodes. Thus, degree(w) — n~ 1 and tt^ — 1/n for all v G V, so we have by Theorem 
[T]that E[T{i)] = t/n. This result should be intuitive by symmetry. Lemma[T]also gives us an upper bound 
on the infection probability, E[7(t)] < 1 — cxp{~(3t/n). Roughly speaking, this says that it takes time of 
order n/p for the susceptible individual to become infected; for t -^n/ (3, the probability of being infected 
is vanishingly small. Again, this is consistent with intuition. 

Regular graphs A graph G = (V, E) is said to be r-regular if degree(u) = r for all v G V . Thus, a complete 
graph is regular with r = rt — 1 . It is readily verified that tTu = 1/n for all w G y if G is r-regular for any 
r > 2. (If r = 1, then G is a matching and is not connected.) Hence, if G is connected, we have the same 
estimates for r(t) and 7(t) as for the complete graph, which is a special case corresponding to r = n — 1. 

The next examples we consider will be families of random graphs widely used in practice to model 
networks. 

Erdos-Renyi random graphs The Erdos-Renyi graph G{n,p) is defined as a random graph on n nodes, 
wherein each edge is present with probability p, independent of all other edges. We consider a family 
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of such random graphs indexed by n, and take p to be a function of n chosen so that np > clogn for 
some constant c > 1. We also condition on the graph being connected. For p as above, the probabihty 
of connectivity tends to 1 as n tends to infinity, so conditioning on connectivity does not aher any of 
the estimates we shah derive later for the coincidence time on such graphs. In this model, the node 
degrees are identically distributed Binomial random variables with parameters {n — l,p). In particular, 
they concentrate around the mean value of np, and have exponentially decaying tails away from this value. 
Thus, while Erdos-Renyi graphs are not exactly regular, they exhibit considerable homogeneity in node 
degrees. 

Power law random graphs In contrast to the above graph models, many real-world networks exhibit 
considerable heterogeneity in node degrees, and have empirical degree distributions whose tails decay 
polynomially; see, e.g., [TJIZ]- This observation has led to the development of generative models for graphs 
with power-law tails [U [2] as well as random-graph models possessing this property [1] . For definiteness, 
we work with the model proposed in |4j , but we believe that similar results will hold for the other models 
as well. 

In the model of [4], each node v is associated with a positive weight Wy, and edges are present inde- 
pendently with probabilities related to the weights by 

P((u, v)eE)^ where W=J2^^- (4) 

xev 

We assume that W ^ ^max' that the above defines a probability. It can be verified that IE[degree(v)] — 
Wy and so this model is also referred to as the expected degree model. The model allows self-loops. 
The Erdos-Renyi graph G{n,p) is a special case corresponding to the choice = np for all v € V. If 
the weights are chosen to have a power-law distribution, then so will the node degrees. The following 
3-parameter model for the ordered weight sequence is proposed in [4J, parametrised by the mean degree 
d, the maximum degree m, and the exponent 7 > 2 of the weight distribution: 

Wi — m(l + — j '' , i — 0,l,...,n—l, (5) 



where 



Note that W = J27=o ^ 

We consider a sequence of such graphs indexed by n. The maximum expected degree m and the 
average expected degree d may, and indeed typically will, depend on n. In models of real networks, we 
can typically expect d to remain bounded or to grow slowly with n, say logarithmically, while m grows 
more quickly, say as some fractional power of n. In this paper, we only assume the following: 

d > S > 0, d = o(m), m < V nd, — — oyni-^ j . (7) 

Here, 5 is a constant that does not depend on n. In other words, the average expected degree is uniformly 
bounded away from zero. The third assumption simply restates the requirement that Wq < W, so that ^ 
defines valid probabilities. The last assumption ensures that iq, defined in ([6|), tends to infinity. 
We now describe our results about these models. 

Theorem 2. Consider a sequence of graphs G = {y,E) indexed by n = \V\. On each graph, consider 
two independent random walks with initial condition Xq, Yq chosen independently from the invariant 
distribution ir for the random walk on that graph. 

We have E[r(i)] = t/n for regular graphs, including the complete graph, on n nodes. 

For Erdos-Renyi random graphs G{n,p) conditioned to be connected, and having np > clogn for some 
c > 1, we have E[r(t)] ^ t/n, as n tends to infinity. 
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Finally, consider a sequence of power law random graphs defined via ^ and and satisfying the 
assumptions in Qj. Then, we have the following: 



nE[r(t)] 



c, if7>^, 
c(logm), i/7 = 3, 
c{m.df~^, i/2<7<3, 



where c > is a constant that may depend on 7, but not on n, m or d. 

3 Proof of Theorem [2] 

If the graph G is regular, then, by ([3]), 7r„ = \/n for all v € V . Hence, the claim of the theorem follows 
from Theorem [T] 

In order to estimate E[r(t)], we need to compute 

V- 2 E.,ey degree(^)^ 



vev (Et,ev degree(t;)) 



Define 



where A^v — v) G E), and 



D=Y^ degree(w) = ^ (9) 



= degree(t;)(degree(?;) - 1) = '^A^iA^ 



'1.^3 

vev 

We will derive the first and second moments of the variables D and D2- It then suffices to use 
Chebyshev's inequality to establish concentration results for both variables D and D2. By ([5]) and Theorem 
[U and the fact that J^vev degree(?;)^ = D2 + D and D = J^vev degree(u), we will have an estimate of the 
coincidence time that holds whp. 

We begin by computing the mean and variance of D in the expected degree model with arbitrary 
weight sequence {wi, z = 0, . . . , n — 1}. 

For notational convenience, we define 



-j^n — 1 ^ n — 1 

W ~ — Wi, w'"' = — ufi , k — 2,3, 



n '■ — ' n 

i=0 i=0 



We obtain Erdos-Renyi graphs G{n,p) by setting Wi = np for all i, and so, w'^ = [np)^ for such graphs. 

Next, consider power-law graphs with weight sequence specified by ([5]) and Since iq tends to 

infinity by assumption, we have for such graphs that 



h. n— 1 

TO 



n ^ — ' V In 

1=0 ^ 



ax 



n Jo ^ ^o 

to'^- / {l + x)- — dx. (11) 
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Now, straightforward calculations yield that w ^ d for all 7 > 2, whereas, for fc > 2, we have 

(^3T)fe?J^rf^ if7>fc+i, 

^l^^dMogf, if7 = fc + l, (12) 

u.^ir-w-.) ^'^"^^'''^^-^ if2<7<fc + i- 

We can now compute the mean and variance of D, the sum of node degrees. 

Lemma 1. Consider a random graph G ~ (V, E) specified by the expected degree model with an arbitrary 
weight sequence {w^^v G V} satisfying W > w^^^, where W — J2vgv^'"- '^^"^ of node degrees, 

D, be defined as in (0). Then, we have 

E[D] = nw, 

VarZ? = 2(nw-[ — )]-(—-— ), 13 

where n = \V\ is the total number of nodes. 

In particular, if G is the Erdos-Renyi random graph G{n,p), then 

E[D] = n^p 

Var(D) (2n-l)np{l-p) ^2n'^p{l-p), (14) 

whereas, if G is a power law random graph satisfying the assumptions of Theorem then 

M[D] = nd while Var(Z?) ^ 2nd, whp. 
Proof. It is immediate from © that 

E[D] = F{{u,v) e E) 



~W 



which estabhshes the first equahty in (fT51) . Next, rewrite © as 

n n n 

2 — 1 j — 2—1 

and observe from the independence of the edges that 

71 n n 

Var(i^) - 4^ ^ Var(Ay) + ^Var(v4,,) 

11 n 11 

= Var(^«j ) - J2 Var(^»»)- 

i—1 j—1 i—1 

Now, Var(A„„) = P{{u,v) G E){l-P{{u,v) e E)), and so, 

v-(^)^^i:n^-1#)-i:(#-^ 

i=l j=l i=l 

Upon simplifying, this yields the second equality in (|13p . Now, using the fact that w'' = {np)^ for 
Erdos-Renyi graphs G{n,p), we readily obtain (fT4|) . 
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Next, suppose G is a power-law graph (more precisely, Gn is a sequence of power law graphs) satisfying 
the assumptions of Theorem [21 It foUows from ([TI]) that 

id' log 1, if 7 = 5, (15) 

(2-2)2~i 



-(P-^ra^--^, if2<7<5. 



while 



(7-1)^-^(5-7) 



(7-l)(7-3)" ' II 7 > «>) 

id^logf, if 7 = 3, (16) 

(,%^r^L,) <P-'m^-\ if2<7<3. 



and w ~ d for all 7 > 2. 

By (fT3|) . it suffices to show that 



{ ] = oind) and ^ — ohi^d) (17) 



in order to show that Var(Z?) ^ 2nw ^ 2nd. 

Suppose first that 7 > 3. Then, by (fTB]) and the fact that W = nd, 

where the last equality follows by and the fac that d < n. 
Now let 7 = 3. Then, by (fTH]) and the fact that W — nd, 



^(^)'-o(^log^)^o(™^log^)=o(l), 
nd \ w / \n a/ \ n m a/ 

where the last equality follows by On the other hand, if 2 < 7 < 3, then, by (fT6)) . 

1 /li^\2 _ ^/d27-5^6-27 

nd \ w J V n 

^ «((^)n-(')> 

where we have used the inequality m < ^/nd from ([7]) to obtain the second equality. To obtain the last 
equality, note that it follows from ^ that m = o{n) and hence that d = o{n) as well. We have thus 
established the first equality in (fTT]) for all 7 > 2. The proof of the second equality is similar and is 
omitted. This completes the proof of the lemma. □ 
The following corollary is now an easy consequence of Chebyshev's inequality. 

Corollary 1. If Gn, n £ N is a sequence either of Erdos-Renyi random graphs or of power-laws random 
graphs satisfying the assumptions of Theorem\^ then the sum of node degrees D concentrates at its expected 
value in the sense that D ^ E[£'] with high probability (whp). 

We now establish a similar concentration result for the sum of squared degrees. To this end, recall that 

Xy = degree(w)(degree(?;) - 1) = AyjAyj 

vev 

We have the following: 
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Lemma 2. Let D2 be defined as in UU\). We then have 

\ w / 



Proof. We first note that 



Var(L>2) < 4nw3 + 2nui2+4n-^^. (18) 

w 



1^2 



up" 



nuP' 



Therefore, 



which is the first part of (fT8)) . Next, for distinct nodes u,v €V , we have 
Cov(X„,X„) = ^^Cov(A„Aj„,^fe„A,^) 

= 4E[A„,„](1 - E[^„,„]) ^ E[A„] ^ E[A,„]. 

The second equality above holds because, by the independence of edges, the indicator random variables 
AiuAju and AkvAiy corresponding to the open triangles (or 2-stars) iuj and kvl are independent unless 
two of the edges are the same; the only way this can happen is if (m, v) is a common edge and there are 
4 possible node labeUings corresponding to each such edge set. Now, recall that E[A„,i,] = WuWv/W and 
X]iIE[Ai„] = E[degree(w)] — Wu- Hence, we see from the above that 

2 2 

0<Cov(X,,X,)<4^. (19) 

Similarly, we obtain 

Var(X„) = YYGoY{Ai^Aju,AkuAiu) 

i^j k^i 

= 4 ^ ^ ^ Cov{A^uAju,AjuAiu) 
+2^Var(A„Aj„) 

+2^E[A„Aj„]. 

Using the fact that distinct edges are independent, we get 

Var(X„) < 4wl + 2wl. (20) 
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Now, by (Uni), (Ull) and (1201) 



w 



Computing the above sums yields the second part of u8\i . □ 
We now specialise the results to Erdos-Renyi and power law random graphs, showing that D2 concen- 
trates near its expected value with high probability. 

Lemma 3. Suppose G{n,p) is a sequence of Erdos-Renyi random graphs indexed by n (where p depends 
on n but this is not made explicit in the notation), and that np is uniformly bounded away from zero. Then 
D2 ^ E[D2] ^ n^P^ whp. 

Proof. We have, by Lemma 12] and the fact that = [np)^ for the Erdos-Renyi random graph G{n,p), 
that 

E[D2] ^ n^(n ~ \)p^ - nV, Var(i:>2) < 8nV + 2nV- 
Hence, by Chebyshev's bound, wc obtain for all e > that, 

^i\D2-nm>enm < 



< 



sm[D2Y 

1 

e^{n — l)^p 
1 



e'^n{n — Vf^p^ 

Now, by the assumption that np is bounded away from zero, {n — and nin — Vf^p^ tend to infinity as 
n tends to infinity. Thus, P(|I?2 — ^[D2\\ > £^[02]) tends to zero for all e > 0. This establishes the claim 
of the lemma. □ 

Lemma 4. Suppose Gn,n G N is a sequence of random graphs satisfying the assumptions in Theorem\^ 
with 7 > 2. Then, D2 ~ IE[-D2] whp, and 

{cnd^, i/ 7 > 3, 

cnd2 (log to), 1/7 = 3, 
cnd'f-^m^-'', i/2<7<3. 

Proof. We will show that Var(Z)2) ~ o(E[Z)2]^), so that the claim follows by Chebyshev's bound, as in 
the proof of the previous lemma. We will consider separately the parameter ranges 7>4, 3<7<4 and 
2 < 7 < 3, where 7 is the exponent in the power law describing the degree distribution. 

In the following, ci, C2, . . . will denote generic positive constants, not necessarily the same from line to 
line. Recall that w ^ d. 

Suppose first that 7 > 4. Then, by ^2]) . — 0{d^ log ^) and ^ cid^. Therefore, by Lemma 121 

E[D2] ^ cmd^ 

and 

Var(i:i2) = o{nd^ log ^ + nd^'^ = o(jid^ 7) ' 

where the last equality holds because of the assumption in ^ that d > S for some constant S > 0. Thus, 
we see that 

Var(L>2) ^/ 1 , m\ 

since m < n. 
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Suppose next that 3 < 7 < 4. Then, by = 0{d''-^m^-^), while w'^ ~ ci(P if 3 < 7 < 4 and 

~ C2(f log ^ if 7 = 3. Therefore, by Lemma [21 

E[D2] > cind? - C2d^ log^ ^ 

> cin(f - C2(flog^n ^ n{n(f), (21) 



whereas, 

Varpa) < cind'^-^m'*-'^ + C2nd2 log ^ + cgndMog^ ^ 
< cirid'''~^m^~''' + C2nd'^ log^ 

a 

We have used the assumption that d is uniformly bounded away from zero to obtain the second inequality 
above. Since we also assumed in ([7]) that d — o{m), we have (d/m)^^'^ log^(r7i/(i) = o(l) for all 7 < 4. 
Hence, Var(Z32) ~ 0{nd^~'^m^~'^). Combining this with ((2T|) . we get 



Var(i?2) ^ ^ 



2J 



We have used ([7]) to obtain the second equality above and the fact that 7 > 3 to obtain the last equality. 
Moreover, E[Z?2] ^ cnd^ for 3 < 7 < 4, whereas E[Z)2] ~ cnd^ log(m) for 7 = 3. 

Finally, suppose that 2 < 7 < 3. Then, by ([12]), w'^ = 0{d''~^m'^~'^) and w'^ ~ cid^~^ni^~'' , so that, 
by Lemma [2j 

Now, by O, {m/df~^ = o(n(3-T)/(')'-i)) = o(n) since 7 > 2. Consequently, 



On the other hand, 



Var(i:>2) < Cind'<-'^m^-'< + C2nd''-'^7rC'~'^ 

^ J7-1 4-7/1 ,C2 

< cina ' m ' I 1 H 

V TO 



C3 - 
TO Vm 



Hence, 



Var(D2) _ q( nd-^'^rn^'^ 



¥.[D2Y \n^d^'>-'^m<^-^-' 



nm \ d 

Now, by ([7|), and the fact that 7 > 2 we have (m/d)''^^ =^ o{n). Since the maximum degree m is assumed 
to grow as a power of n, we have ^^^^3^2"' — o{\). Note that E[D2] ~ cnd''~'^Tin?~'^ , for 2 < 7 < 3. 

Using Chebyshev's inequality, this establishes the claim of the lemma. □ 
To complete the proof of Theorem [2l it suffices to combine the results of lemma [T] and lemma [4| with 
the fact that Y.v '^v = 
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4 Conclusion and future work 



In this work we have presented a simple model for the spread of epidemics where individuals are mobile. 
In this framework we were interested in the setting where there are two individuals one infected and one 
healthy both performing random walks on the network. Our preliminary investigation highlights the effect 
of the topology on the spread of an epidemic, motivated by networking phenomena such as worms and 
viruses, failures, and dissemination of information. Under this natural model, we provided an explicit 
relationship between the structure over which the walks are performed and the coincidence time of the 
two walkers. To this end we analysed both homogeneous (regular, complete and Erdos-Renyi graphs) and 
heterogeneous (power-law graphs) networks. We pinpointed the existence of a phase transition for the 
coincidence time in the case of power-law networks depending on the parameter of the power-law degree 
distribution. We also derived bounds on the probability of infection. 

As a final remark, we propose some several interesting directions to pursue the work presented here. In 
our present model individuals are supposed to start their walks in stationary regime. This can be relaxed 
since the networks we study are expanders and thus random walks on such networks have nice mixing 
properties as illustrated in [S] through the computation of the isoperimetric constant of the underlying 
graphs. We also anticipate that similar results can be derived when considering k walkers as long as k is 
small with respect to n the number of sites in the network. 
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